Collective sedimentation of colloids in grafted polymer layers 
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We present Monte Carlo simulations of colloidal particles pulled into grafted polymer layers by 
external fields. The insertion free energy of a single colloid into the polymer layer is qualitatively 
different for surfaces with an ordered and a disordered distribution of grafting points. Moreover, 
the tendency of colloidal particles to traverse the grafting layer is strongly size dependent. In dense 
colloidal suspensions, under the influence of sufficiently strong external force, a collective instability 
allows the colloids to penetrate and form internally ordered, columnar structures spanning the 
polymer layer. Depending on the conditions, these colloidal clusters may be isolated or laterally 
percolating. The morphology of the observed patterns can be controlled by the external fields, which 
opens up new routes for the design of thin structured films. 



Grafted polymer layers can prevent the deposition of 
colloidal particles on solid surfaces, which is exploited 
in various applications like colloidal stabilisation, anti- 
fouling surfaces [T| and in biomedicine [21 [3] . Biological 
surfaces such as the lining of the intestine or the blood 
vessels walls are coated with polymers that inhibit the 
absorption of too large particles. Moreover, membrane- 
attached polymers over-expressed in cancer cells present 
a barrier for conventional drug delivery processes [1]. In 
several other applications, the insertion of colloids or 
nano-particles into polymer layers is exploited in order to 
fabricate functional responsive materials [5]: plasmonic 
nanostructures [6j convert chemo-physical forces into op- 
tical signals; ultrafast switches [7], organic memory de- 
vices and advanced photovoltaics |9l [10] are based on 
incorporating metallic particles into polymers. 

The physical properties of grafted polymers layers are 
governed by the configurational entropy of tethered poly- 
mers and have been studied extensively for different 
grafting and solvent conditions [TTlfM] . The penetra- 
tion of particles into polymer layers has been recently 
addressed in theoretical and numerical studies [151 [H] . 
The insertion free energy of an isolated protein [T71 [TS] 
or colloid [El [20] has been explored for varying graft- 
ing densitiies, particle sizes and shapes [H]. Colloids 
soluble in polymers penetrate up to a depth determined 
by their size [2?| . while insoluble particles form aggre- 
gates, which are expelled from the brush after reaching 
a critical size [23]. The polymer-mediated colloidal in- 
teractions are typically of the order of fc^T [2T1[53] (i.e., 
comparable to the entropic terms). The colloid-polymer 
mixtures are thus inherently less ordered than molecu- 
lar systems. External forces could potentially be used to 
re-gain the order: due to the constrained translational 
degrees of freedom of the colloids pulled into the layers, 
the weak polymer-induced interactions could steer the 
assembly into ordered patterns. However, little is known 
about the collective particle sedimentation in a brush un- 
der the influence of a constant force. 

Here we report Monte Carlo simulations of polymer- 
insoluble colloidal particles in grafted polymer layers. 
The system comprises of polymer chains end-grafted to 
a flat surface and hard sphere-colloids with diameter a 
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FIG. 1. Isolated colloids in polymer layers with grafting den- 
sity p = 1.0 and number of blobs per chain Ip — 20. (a) 
The insertion free energy profile Fp{z) for three realizations 
of the surface grafting: quenched disorder (black solid line), 
square crystaline order (red solid line), and annealed disorder 
(green dashed line), (b) The distribution of (2D) hole sizes for 
the three scenario. The solid lines represent the analytically 
derived expressions (see Online SI), symbols are numerically 
determined. 



immersed in a neutral good solvent. The colloids inter- 
act via hard-core and polymer-mediated interactions and 
are additionally subject to a constant external force fg, 
whose significance is measured by the ratio of potential 
versus thermal energy: g' = fga/ksT. A possible exam- 
ple of the external force is gravity: fg ~ g{pc — Ps)7rcr^/6, 
with pc — Ps the excess colloidal density in the solvent 
and g the standard gravitational acceleration. For grav- 
ity, the ratio g' oc a'^ strongly depends on the colloidal 
size: it is efficient for large-enough (micron-size) colloids 
and does not affect nano-particles or polymers in solu- 
tion. Polymers can thus be treated as self-avoiding ran- 
dom walks. We follow a coarse-grained model [25^ where 
the self-avoiding polymer chains are represented by Ip 
soft blobs with radius of gyration rf, [26j. The grafting 
density at the surface, p = Np/R^, determines whether 
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the layer is in the dilute "mushroom" (p < 1) or in the 
dense "brush" scaling regime (p > 3) [271 [28]. We studied 
the transition regime with the mean spacing between the 
grafting points roughly similar to the polymer Rg |29| . 
Expressing the grafting density as number of anchors per 
colloidal diameter squared: p = i5B^/a^, we focused on 
the regime 0.5 < p < 5, where the probability of inserting 
a colloid with a ^ Rg is non- vanishing. 

The free energy profile F{z) of a single colloid is the 
sum of its potential energy in the external field and the 
insertion free energy Fp(z): 

/3F{z)=g'^+l3Fp{z) . (1) 

The potential energy term is trivial and decouples from 
the polymeric degrees of freedom, therefore the insertion 
free energy is the key quantity governing the penetra- 
tion of colloids into the polymer layers. We have used 
the above described coarse-grained model of polymers 
and evaluated the insertion free energy Fp{z) by the 
Wang-Landau technique [30] for a wide range of poly- 
mer grafting densities p and number of blobs per chain 
Ip. In Figure 1(a) we compare polymer brushes char- 
acterized by ordered and disordered arrangement of the 
grafting points: Fp{z) is plotted for layers with identical 
p but different spatial distribution of the grafting points: 
"quenched disorder" (random Poissonian process), "an- 
nealed disorder" (random grafting points relaxed before 
the simulation), and "order" (grafting points on a square 
lattice). The insertion free energy is qualitatively dif- 
ferent in ordered and disordered case: the repulsion is 
strongest in case of ordered grafting where the free energy 
is monotonically decreasing with the height. In disor- 
dered case the free energy features a barrier. The height 
and position of the barrier and the depth of the mini- 
mum at the surface can be controlled (see Online SI), 
which enables reversible adsorption and slow release of 
the colloids at disordered surfaces. 

The emergence of the free energy barrier (also ob- 
served in other simulations [151 I18| ) is surprising: the 
mean monomer concentration profile is equal for ordered 
and disordered grafting and, at p = 1.0, it is monotoni- 
cally decreasing with height (see Online SI and [HI HI]). 
The insertion free energy is obtained by exponential av- 
eraging of such profiles and is determined by both, the 
mean monomer concentration and its fluctuations. The 
fluctuations are characterized by the probability distri- 
bution of holes (regions void of monomers) [33) . In the 
inset (Figure 1(b)) the 2D hole-size distributions at the 
bottom surface are plotted: in case of disordered graft- 
ing characterized by the Poissonian statistics, there is 
an appreciable probability of finding large holes, while 
the ordered surfaces are characterized by small density 
fluctuations with a cut-off in the hole size. Close to the 
grafting surface the monomer fluctuations are "frozen" 
because the anchoring points are immobile, while away 
from it they relax: the surface-imposed effects vanish 
at a height equal to the characteristic length-scale [53]. 




FIG. 2. (a) Total free energy ([T| of sedimenting colloidal 
particles obtained by simulations at Ip = 40 and p — 1.0 
for silica colloids with three sizes: ai — 0.67/im (red line), 
a2 ~ 1.00/im (black line), and o"3 — 1.33pm (green line). 
Respectively, g2 — 2.0 (corresponds to micrometer silica par- 
ticles in water), g'l — 0.40 and g'g — 6.3. The dashed lines 
represent the insertion free energy Fp{z). (b) Height distri- 
bution P{z) oc e"''^'^-' for the same system. The histograms 
illustrate the probability to find the colloids adsorbed to the 
surface or above the brush (left set of histograms: z < Rg 
and right set: z > Rg)- (c) Snapshots from Monte Carlo 
simulations of dilute colloidal suspensions at conditions cor- 
responding to (a) and (b). 



The free energy cost of inserting a particle into a dis- 
ordered brush therefore initially grows as a function of 
height and eventually decreases along with the average 
monomer concentration - resulting in the barrier with a 
maximum around z Rg. The barrier height scales with 
p and with ct^; in our system its typical magnitude is 
several 10 ksT (see Online ST). 

The potential energy term ([T]) significantly alters the 
free energy landscape: for large enough g' the unstable lo- 
cal minimum becomes stable and the colloids are strongly 
confined to the surface. Since in the case of gravity g' sen- 
sitively depends on the colloidal size, this property can be 
used for particle-size sorting. In Figure 2 we plot the total 
free energy profiles (a) for colloids of three different sizes 
together with their height probability distribution (b). 
Counterintuitively, the smallest colloids gather on top of 
the brush (the gravity is too week to overcome the steric 
barrier) , while slightly larger colloids are distributed in a 
bimodal way and even larger ones are strongly confined to 
the bottom surface. This single-colloid picture also gov- 
erns the density profiles in semi-dilute systems as seen in 
the simulation snapshots in Figure 2c)). The amount of 
particles drawn into the polymer layer reflects the bal- 
ance between the polymer-induced pressure pushing the 
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particles away and the external forces pulling them to- 
wards the surface. At a critical grafting density p — po 
the two effects are balanced: for p > p^ no particles 
can penetrate the polymer layer [34j . For p < pQ the 
particles are pulled in the polymer layer by the excess 
external pressure - thereby displacing the polymers and 
increasing their effective density. When many colloids 
are inserted, the polymers mediate effective many-body 
interactions among the particles and the insertion free 
energy is no longer described by the single-particle term 
Fp{z). In general it depends on the positions of all parti- 
cles. To explore the multi-particle systems we have per- 
formed Grand-canonical Monte Carlo simulations with 
colloidal particles coupled to a reservoir with a fixed 
chemical potential at various polymer grafting densities, 
chain lengths, particle concentrations and external forces. 

At large-enough pulling g' the penetrating particles 
start forming crystalline clusters at the bottom surface 
that grow towards the top of the brush [35|. Figure 3 
shows typical growth scenarios on disordered and ordered 
surfaces. In the initial stages there is a difference: the 
shape of the clusters above the ordered surfaces is nar- 
rower at the bottom (resembling a table-top) due to the 
strong repulsion of colloids from the surface, while in case 
of disordered grafting the local free energy minimum en- 
ables pyramid-like structures. In both cases, however, 
similar final structures emerge: clusters with a uniform 
vertical profile spanning the brush from bottom to top 
resembling straight cylindrical towers. The maximum 
load of colloids in the polymer layer is achieved when the 
effective polymer density in the space around the col- 
loids equals Po- Consequently, the thickness of a fully 
loaded layer cx Pq^^ [TT] does not depend on the poly- 
mer grafting density and can be controlled by inserting 
particles into the polymer layer (experimentally, an ex- 
cessive amount of colloids could be added and the ones 
above the brush could be simply washed-off). Given the 
uniform vertical profiles observed in the colloidal clusters 
(Figure 3), the particle load is directly correlated to the 
colloidal surface coverage 77 of the two-dimensional hor- 
izontal cross-sections (top- view snapshots in Figure 4): 
77 = Scoi/S = I — where Scoi/S is the fraction of the 
area on the snapshot covered by the (red) colloids. 

In Figure 4 we present the phase diagram as a function 
of the grafting density p and the critical density pq oc g' . 
We have investigated both, ordered and disordered sur- 
faces and observed no important differences in the lateral 
morphology. Different structures observed in the simula- 
tions are depicted by symbols and their phase boundaries 
are approximately marked by the solid lines, which are 
lines of constant colloidal surface coverage rj. The white 
lower region of the diagram is the "super-critical" graft- 
ing region where no colloids penetrate the brush and no 
ordered structures are formed. The grey-shaded region 
depicts the regime where the colloids can penetrate the 
brush but the effective polymer density is too small and 
the entropy prevents particle ordering [36] . In order to 
understand the micro-phase separation leading to these 
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FIG. 3. Cluster growth as a function of the number of colloids 
in the system. Sedimentation sequence on ordered (left-hand 
side) and disordered (right-hand side) surfaces is depicted. 
Initially, "pyramide" and "table-top" structures are formed 
on disordered and ordered surfaces, respectively. Eventually, 
in both cases "critical" clusters with uniform vertical profiles 
spanning the polymer layer are formed. The number of blobs 
per chain is Ip — 40, the anchoring density p = 4.0, and the 
effective gravity g' — 6.0. 



various patterns we consider two mechanisms: the deple- 
tion of polymers around solid objects that favours macro- 
phase separation, and an elastic-like penalty due to the 
fact that the polymers are grafted that promotes colloid- 
polymer mixing. The competition between such terms, 
combined with the action of the external forces, results 
in the formation of finite and ordered lateral patterns 
uniquely characterized by the colloiodal surface coverage: 
at supercritical conditions (77 = 0) the colloids are in a 
liquid state on top of the brush; at small rj isolated circu- 
lar clusters emerge. Upon increasing the colloidal surface 
coverage snake-like objects, percolating structures and 
finally inverted clusters are formed. The morphologies 
are quite robust and are characterised by a well-defined 
length-scale (the width of the towers, snakes, or percolat- 
ing structures). Interestingly, this length-scale does not 
depend on the grafting density but scales with Ip (i.e. 
with the polymer size) and is roughly equal to the brush 
height. The stability and properties of the reported pat- 
terns is reproduced by a phenomenological theory pre- 
sented in Online SI. 

Our results provide an insight into many-body 
polymer-induced particle interactions, which are of key 
relevance for understanding various biological processes, 
and to design novel materials and biomedical applica- 
tions. Metal-coated micro-particles pulled into grafted 
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FIG. 4. Phase behaviour as a function of the polymer graft- 
ing density and the external pressure. Various lateral patterns 
are depicted by symbols: isolated cylindrical towers (blue cir- 
cles), "snakes" (blue triangles), percolating structures (red di- 
amonds) and inverted towers (red squares). The black lines of 
constant surface coverage rj separating the phases are guides 
to the eye. Top-view snapshots from the simulations display- 
ing the lateral morphology are shown at the side and marked 
by the same symbols. The snapshots are taken at an interme- 
diate height within the polymer layer and correspond to the 
value of the critical density po = 6.0 with g' = 6.0, Ip — 40 
and the grafting density ranging from p = 1.0 (inverted tow- 
ers) to p = 5.0 (isolated towers). The typical length-scale and 
the strong crystalline order are clearly revealed. 



layers by gravity can form thin metallic films with unique 
mechanical, electric conductivity or heat transport prop- 
erties. If mechanical stress can be applied to stretch or 
bend the substrate, the grafting density and the mor- 
phology of the emerging structures can be varied in a 
controlled way, which should provide a useful insight for 
surface characterization. We furthermore predict that 
grafted surfaces can be used as a sorting device for poly- 
disperse particles or particle aggregates, as well as to con- 
trol the rates of chemical reactions. It has been shown 
[T5] that the polymer brushes can affect protein folding 
processes, which is not surprising in light of our results. 
In order to fabricate the predicted multi-particle struc- 
tures with nano- instead of micro-particles, ultra cen- 
trifuge can be used instead of gravity. The morphology 
of the reported structures is strongly desirable in mod- 
ern nano-responsive applications like photovoltaics [37], 
smart glasses [7] and miniature sensors [38] . 
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